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ABSTRACT 

Earth-mass "microhalos" may be the first objects to virialize in the early universe. Their ability to 
survive the hierarchical clustering process as substructure in the larger halos that form subsequently 
has implications for dark matter detection experiments. We present a large N-body simulation of 
early substructure in a supersymmetric cold dark matter (SUSY-CDM) scenario characterized by an 
exponential cutoff in the power spectrum at M c = 1CP 6 M©. The simulation resolves a 0.014 M Q 
parent "SUSY" halo at z = 75 with 14 million particles. On these scales the effective index of the 
power spectrum approaches —3, and a range of mass scales collapse almost simultaneously. Compared 
to a z = galaxy cluster substructure within our SUSY host is less evident both in phase-space and 
in physical space, and it is less resistant against tidal disruption. As the universe expands by a 
factor of 1.3, we find that between 20 and 40 percent of well-resolved SUSY substructure is destroyed, 
compared to only ~ 1 percent in the low-redshift cluster. Nevertheless SUSY substructure is just 
as abundant as in z = galaxy clusters, i.e. the normalized mass and circular velocity functions 
are very similar. The dark matter self-annihilation 7-ray luminosity from bound subhalos and other 
deviations from a smooth spherical configuration is at least comparable to the spherically-averaged 
signal in the SUSY host, and at least three times larger than the spherically-averaged signal in the 
cluster host. Such components must be taken into account when estimating the total cosmological 
extragalactic 7-ray annihilation background. The relative contribution of bound substructure alone 
to the total annihilation luminosity is about four times smaller in the SUSY host than in the z = 
cluster because of the smaller density contrast of sub-microhalos. 

Subject headings: methods: N-body simulations - methods: numerical - early universe - structure 
formation 



1. INTRODUCTION 

The key idea of the standard cosmological paradigm 
for the formation of structure in the universe, that pri- 
mordial density fluctuations grow by gravitational insta- 
bility driven by cold, collisionless dark matter (CDM), is 
constantly being elaborated upon and explored in detail 
through supercomputer simulations, and tested against a 
variety of astrophysical observations. The leading candi- 
date for dark matter (DM) is the neutralino, a weakly in- 
teracting massive particle predicted by the supersymme- 
try (SUSY) theory of particle physics. While in a SUSY- 
CDM scenario the mass of bound DM systems (halos) 
may span about twenty order of magnitudes, from the 
most massive galaxy clusters down to Ea rth-mass clumps 
fe.g. lCreen. Hofmann. fc Schwarzll2004l) . it is the small- 
est DM microhalos that collapse first, and only these 
smallest scales are affected by the nature of the relic DM 
candidate. 

Recent numerical simulations of the collapse of the ear- 
liest and smallest gravitationally bound CDM clumps 
(piemand. Moore, fe Stadelll2005l iGao et alJl200H) have 
shown that tiny virializcd microhalos form at redshifts 
above 50 with internal density profiles that are quite 
similar to those of present-day galaxy clusters. At these 
epochs a significant fraction of neutralino DM has al- 
ready been assembled into non-linear Earth-mass over- 
densities. If this first generation of dark objects were 
to survive gravitational disruption during the early hi- 
erarchical merger and accretion process - as well as 
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late tidal disruption from stellar encounte rs ( Zha o et alJ 
l2005tlBerezinskv. Dokuchaev. fc Eroshenkol200fifl - then 
over 10 15 such clumps may populate the halo of the Milky 
Way ijDiemand et all2005[) . The nearest microhalos may 
be among the brightest sources of 7-rays from neutralino 
annihilation. As the annihilation rate increases quadrat- 
ically with the DM density, small-scale dumpiness may 
enhance the total 7-ray flux from nearby extragalac- 
tic systems (like M31), making them detectable by the 
forthcoming GLAST satellit e or the next-generation of 
air C erenkov telescopes (e.g. lOda. Totani. fe Nagashimal 
l2005|) . Dark matter caustics might occur within the first 
generation of halos. Caustics are not resolved in N-body 
simulations like the ones discussed here and might have 
significant effects on DM d etection experiments (e.g. 
IMohavaee fc Shanda,rmll200l. 

The possibility of observing the fingerprints of the 
smallest-scale structure of CDM in direct and indi- 
rect DM searches hinges on the ability of microhalos 
to survive the hierarchical clustering process as sub- 
structure within the larger halos that form at later 
times. In recent years high-resolution N-body simula- 
tions have enabled the study of gravitationally-bound 
subhalos with M^hh/ Mhain >. 10~ 6 on galaxy (and 



galax y cluster) scal e s llMoore et al.l 119991 iKlvpin et al 
IStT^nTetTaTI )200lFlPiernai 1 d. Moore, fc Stade 
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~IOa,o et a,1.ll2004l IReed et a,l.ll2005tt. The main dif- 
ferences between these subhalos - the surviving cores of 
objects which fell together during the hierarchical assem- 
bly of galaxy-size systems - and the tiny sub-microhalos 
mentioned above is that on the smallest CDM scale the 
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effective index of the linear power spectrum of mass den- 
sity fluctuations is close to —3. In this regime typical 
halo formation times depend only weakly on halo mass, 
the capture of small clumps by larger ones is very rapid, 
and sub-microhalos may be more easily disrupted. Yet 
previous simulations that have resolved sub halos above 
1000 M Q down to z = 4 llMoore et al.l|200lD and above 
10 M Q down to z = 49 l|Gao et al.ll2005|) have found that 
the substructure abundance in early DM halos is quite 
similar to that of their present-day counterparts. 

In this paper we extend this range further down to- 
wards the free-streaming cutoff mass of 10 -6 M Q by 
using a mass resolution of 10~ 9 M Q in a large high- 
resolution TV-body simulation of early substructure in 
SUSY-CDM. The simulation resolves a 0.014 M Q host 
halo at z = 75 with 14 million particles within its virial 
radius. The mass resolution and power s pectrum are the 
same as in one of the runs presented in iDiemand et alJ 
lj2005fl . which focused on a region of average density. 
The new simulation discussed in this work has a high 
resolution region that covers an eight times larger ini- 
tial volume [120 (comoving pc)] 3 and is centered on the 
highcst-er peak of the same 3 comoving kpc parent sim- 
ulation cube. 

The plan of the paper is as follows. In § 2 we describe 
the initial conditions, numerical methods, and parame- 
ters of our simulations. In § 3 we present the properties 
of our SUSY parent halo and its substructure, comparing 
them to a low-redshift massive galaxy cluster simulated 
at comparable resolution, and discuss the relative abun- 
dance and survival probabilities of substructure in both 
simulations. In § 4 we estimate the dark matter self- 
annihilation 7-rays signal from substructure at z = 75 
and z = 0. Finally, we summarize our conclusions in § 5. 
Details about finding early CDM substructures both with 
SKID and a new phase-space friends-of-friends (6DFOF) 
algorithm are given in the Appendix. 

2. N-BODY SIMULATIONS 

The results of this paper are based on a comparison 
between two dark matter simulations on vastly different 
scales. One simulation ("SUSY") resolves DM substruc- 
ture all the way down to the free-streaming scale for 
a generic supersymmctric particle of mass mx = 100 
GeV, and is designed to model the formation of the 
smallest and earliest dark microhalos in the universe. 
The other simulation ("B9") models the formation of a 
5.9 x 10 M Q clust er at the present epoch, and was pr e- 
viously discussed in IDiemand. Moore, fe Stadell l)2004|) . 

Both simulations were performe d with PKDGRAV 
<|Stadell l200lt iWadslev et all l200l and employed the 
m ultiple mass particle gr id initial conditions as described 
in IDiemand et all l)2004|) . The SUSY run has a comov- 
ing box size of 3 kpc and wa s initialized at z = 456 
using the GRAFICS package l|Bertschingerl 120011) with 
a DM power spectrum P{k) for a 100 GeV neutralino 
(kinetic dec o upling occurs at Td = 28 MeV) given in 
IGreen et al.l 1)20041) 1 . This power spectrum is close to 
P(k) ~ k~ 3 and exhibits a free streaming cutoff at 
k = 1.65 (comoving pc) -1 , corresponding to a mass scale 
of 1.1 x 10~ 6 M Q . The initial coupling between DM and 

1 As in IDiemandet aT] OoMl we use model "A" of the first 
preprint version of lGreen et alj 120041) jastro-ph /0309621 tl). 



TABLE 1 
Simulation parameters 



%i %f e ^hircs ^p, hires 

(comoving units) (Mq) 
455.8 75.1 0.00216 pc 6.4 x 10 7 0.98 x KT 9 
40.3 0.0 4.8 kpc 2.8 x 10 7 5.2 x 10 7 

Note. — Initial and final redshifts Zi and Zf, (spline) softening 
length e, and number A^ires ar >d masses m Pi hires of high resolution 
particles, for the SUSY (top) and B9 (bottom) runs. 



radiation produces acoustic oscillations that leads to a 
sli ghtly higher cutoff mass o f 4.5 x 10~ 6 M Q (see eq. 26 
in lLoeb fc Zaldarriagal20 05). but we have neglected this 
correction in the present work. 

The SUSY simulation is centered on a rare high-cr 
peak, which is covered by a 120 comoving pc high- 
resolution region consisting of 64 million DM particles. 
This results in a mass per particle of 9.8 x 10 _10 Mq, 
which ensures that even objects at the free-streaming 
limit are resolved with ~ 1000 particles. At z = 75 
the most massive halo in the box has a virial mass 2 of 
-Mhaio = 0.014 M Q , which corresponds to a 3.5-ct fluctua- 
tion. We have adopted the standard concordance ACDM 
cosmology with parameters ttnj — 0.268, f^A = 0.732, 
Hq = 71 kms _1 Mpc _1 , n = 1, and a power spectrum 
normalized to erg = <j{z = 0, r = 8 hr l Mpc) = 0.9 in the 
SUSY and 0.7 in the B9 run. The parameters of both 
simulations are summarized in Tabic ^ 

3. HALO SUBSTRUCTURE 

The fc~ 3 dependence of the power spectrum near the 
SUSY cutoff translates into a a(M) (the linear theory 
rms density fluctuations extrapolated to the present-day 
in a sphere containing a a mass M), that is nearly in- 
dependent of mass. This implies that structures are si- 
multaneously collapsing over a wide range of masses and 
has important consequences for the nature of DM (sub-) 
structure near the SUSY cutoff, as we show below. 

3.1. Structure formation on SUSY scales 

In linear theory the collapse time of a typical spherical 
top-hat overdensity of mass M is given by o-(M)D(a) = 
5 C , where D(a) is the growth factor and 8 C is the critical 
overdensity for collapse. Assuming a flat universe with 
£1m = 1 at these early epochs, a typical 10~ 6 M Q micro- 
halo will form at z = 30, whereas a much more massive 
0.01 M Q parent will collapse at z = 23, i.e. structure for- 
mation on very different mass scales happens almost si- 
multaneously. By contrast, a(M) is significantly steeper 
on galaxy and cluster scales, which leads to a orderly 
hierarchical buildup of cosmic structure. 

The sequence of nested N-body rcsimulations of 
iGao et aTl i|2005l) already showed that large-scale struc- 
ture at z ~ 50 is qualitatively different from that in 
the low-redshift universe. A comparison between the ex- 
tremely rapid assembly of our SUSY halo at redshifts 
130 < z < 75 and the canonical hierarchical growth of 
cluster B9 at late epochs is given in Figure ^ Here we 

2 The "virial" mass is denned to enclose a mean density 200 
times the mean background density in the high rcdshift SUSY run. 
For the 2 = cluster, it is 368 times the mean background density. 
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Fig. 1. — Local DM density maps. The left-hand panels illustrate the almost simultaneous structure formation in the SUSY run at 
different epochs, within a sphere of physical radius r including a mass of 0.014 Mq (equal to Mh a i at z = 75). The galaxy cluster halo B9 
(right-hand panels) forms in the standard hierarchical fashion: the DM distribution is shown within a sphere of radius r including a mass 
of 5.9 X 10 14 Mq (equal to M[ la i at z = 0). The SUSY and cluster halos have concentration parameter [defined as c = 2.15(r v j r /r Vc max ), 
which is equivalent to c = r vlr /r s for a Navarro-Frcnk- White (1997, hereafter NFW) profile] c = 3.7 and c = 3.5, respectively. In each 
image the logarithmic color scale ranges from 10 to 10 6 p cr it (z). 



show maps of the DM distribution at different epochs 



within a sphere enclosing the final virial mass of the 
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Fig. 2. — Phase-space density map for the z = 75 SUSY halo (left) and the z = galaxy cluster (right). Local phase-space densities are 
estimated using the SMOOTH program (available at www-hpcc. astro. Washington. cdu/tools/smooth. html) over 32 nearest neighbors. Note 
the different color scales: relative to the average virial phase-space density (p/fr 3 D ), the logarithmic color scale ranges from 10 to 10 5 in 
the SUSY halo and from 10 to 10 7 in the cluster halo. 



parent halo, using identical (logarithmic) color scales for 
direct comparison. The relative sizes of the most mas- 
sive progenitors are obviously much smaller for the SUSY 
host than for the cluster. The virial mass of the SUSY 
halo grows by a much larger factor than B9 during the 
comparable fractional increase in scale factor shown here. 
Another consequence of nearly simultaneous structure 
formation is that filaments are very large relative to the 
size of the bigge st virialised systems in the same region 
l|Oa.o et al]l2005h . 

The difference in formation histories discussed above 
should be reflected in the density contrast between sub- 
structure and the host, since the average density at viri- 
alization is proportional to the mean density of the uni- 
verse at the time of collapse. The density contrast scales 
approximately as the inverse cube of the scale factor at 
collapse, (ahaio/asubh) 3 - For early microhalos of masses 
10 -6 M Q and 10 -2 M Q this ratio is only 2.5, much smaller 
than the typical contrast ~ 20 of galaxy-scale substruc- 
ture in a cluster host. The expectation is then that 
substructure in B9 will have significantly higher contrast 
than in the SUSY halo 3 . Visual inspection of Figure ^ 
clearly reveals this effect. 

Both the SUSY run at z = 115 and the B9 run at 
z = 0.8 correspond to epochs when the linear growth fac- 
tor is about 0.66 of the final value at z = 75 and 2 = 0, 
respectively. At z = 115 we find 1931 virialized micro- 
halos, less than the 2906 galaxy-size halos and subhalos 
in the cluster forming region at z = 0.8. The abundance 
of bound halos in the SUSY run becomes similar to B9 
only around z as 100. The fractional size of the most 
massive progenitor is also smaller in the SUSY run: at 

3 This simple argument is strictly valid only for typical la fluc- 
tuations of the density field. Subhalos are likely to have a slightly 
larger contrast, since they form from rarer fluctuations that are 
superposed on their parent halo density peak. 



z = 115 this has a mass of 0.01 Mh a io(z = 75), while the 
largest cluster progenitor at z = 0.8 has already a mass 
of 0.25 M halo (z = 0). It takes until z ks 90 before the 
SUSY main progenitor has also grown to a virial mass 
of about one quarter of the final system. This illustrates 
the compression of host- and subhalo formation times 
resulting from the flatness of a(M) on extremely small 
scales. 

Subhalos have smaller internal velocity dispersions 
than the hot uniform density background. Therefore sub- 
halos stand out much more clearly in phase-space density 
p/oy£) than in real space density, as shown in Figure [21 
The phase-space map reveals that the SUSY halo indeed 
has rich substructure at z = 75, comparable in abun- 
dance to the z = B9 cluster. Yet as in physical density, 
the phase-space density contrast is lower in the SUSY 
run compared to the B9 halo. Note the different color 
scales in Figure [21 the scale spans six decades in the B9 
image, but only four decades in the SUSY image. 

3.2. Substructure mass function 

We have constructed substructure catalogs for the 
SUSY and B9 hosts at various stages in their forma- 
tion. We employ t wo different algorithms to identify 
subhalos: a) SKID (jStadelll200l|) . which gives reliable 
subhalo masses by iteratively removing unbound par- 
ticles, but is computationally expensive; and b) a new 
phase-space friends-of- friends method (6DFOF), which 
only produces halo centers and maximum circular veloc- 
ities, but is significantly faster than SKID and is designed 
to take advantage of the higher substructure contrast in 
phase-space. Both techniques are described in more de- 
tail in Appendix [S] Because of the reduced contrast at 
high redshift, the SUSY sub-microhalo catalogs depend 
strongly on our choice of identification method and pa- 
rameters, much more so than in a z = cluster like B9. 
Nevertheless we find reasonable agreement between the 
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two methods when SKID is run with a linking length of 
0.0216 comoving pc and a smoothing scale of 1024 near- 
est neighbors (cf. Appendix |A"|) . 

We define a dimcnsionless subhalo mass, Af su bh/-^haio, 
and plot in Figure 01 the abundance of SUSY sub- 
microhalos at z = 86 and z = 75 against this normal- 
ized mass. The normalized mass function lies close to 
and within the scatter of the average normalized subhalo 
mass function found in galaxy clusters at the present 
epoch, a self-similarity that has been previously found 
to hold over a wide range of mass scales and redshifts, 
from present-day c lusters to dwarf galaxy halos at z = 4 
ijMoore et al.l2001l), and ev en for 10 5 M Q host minihalos 
at z = 49 ijGao et al.ll2005|) . Our simulations extend this 
range down to 10~ 2 M Q SUSY hosts at z > 75: the sub- 
structure populations of young, low-concentration halos 
are approximate scaled copies of each other, independent 
of parent mass. Note that the same self-similarity also 
appears in the normalized distribution of subhalo maxi- 
mum circular velocities (see Figure 0J. 
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Fig. 3. — Differential SUSY sub-microhalo abundance as a 
function of scaled sub-microhalo mass, M su bh/^halo: at z = 86 
(solid) and at z = 75 (dashed). Halo masses were determined 
by SKID using s = 1024. The dash-dotted line corresponds to 
dN/dln(M auhh /M halo ) = 0.01(M aubh /M halo )-\ which is a good 
appro ximation to the z = galaxy cluster relation IDiemand et al.l 
120041) . The dotted lines show the scatter 0.23 dex) around this 
relation. The first bin begins at 256 bound particles. Bars show 
the Poisson errors. The arrow shows the upper limit for the abun- 
dance of sub-M c clumps (128 to 202 bound particles) at z = 86 
estimated by using SKID with s = 64 (see Appendix for details). 

In Figure |S] we present the non-normalized circular 
velocity function for our set of SUSY sub-microhalos. 
The circular velocity distribution for subhalos of 20 — 



300 kms in hosts of mass 10 



it 



10 14 at redshifts 



from zero to two is fitted by ijReed et al.ll2005|) 

dn 



1.5 x 10 



8 w -4.5 



(h 3 Mpc- 3 km _1 s) , (1) 




Fig. 4. — Differential SUSY sub-microhalo abundance as a func- 
tion of scaled maximum circular velocity, t) C| max/«c max, host , at z = 
86 (solid) and z = 75 (dashed). Halos were identified with a phase- 
space friend-of- friend finder (6DFOF). The dash-dotted line is the 
universal normalized circular velo city function of su bstructure in 
hosts of 1011 to 10 14 M found bv lHeed et"aT] 120051) at < z < 2, 

djV/d('Uc,max/-Uc,max,host) = (1/8) (» c ,max /«c, max, host ) ~ 4 • The re- 
ported halo-to-halo variance in the scatter around the universal 
circular velocity distribution is a factor of 2 to 4, the dotted lines 
illustrate deviations by factors of four. 



tion at z = 86 lies right on the extrapolation (over four 
orders of magnitude) of equation Q . The sub-microhalo 
abundance appears to decrease with cosmic time until the 
last analyzed snapshot at z = 75, where it is still within 
the sc atter around the universal function of iReed et alJ 
l|2005() . Unfortunately we cannot quantify the substruc- 
ture abundance at lower redshifts since the parent halo is 
growing very quickly and starts accreting material from 
the lower resolution regions, in the form of heavier par- 
ticles. Note that the evolution in the subhalo mass or 
circular velocity functions should not be used to con- 
strain the survival probability of substructure, since the 
mass of the host (and hence the volume probed) increases 
significantly with time. 

3.3. Formation of halos below the cutoff scale 

The SUSY-CDM scenario adopted here has an expo- 
nential cutoff in the power spectrum P(k) at a scale cor- 
responding to a mass of M c = 10~ 6 M Q . Below this 
mass the formation of structures is sup pressed relative 
to a C DM scenario without such a cutoff (|Diemand et alJ 
120051) . But the suppression is not total: some halos and 
subhalos with masses well below M c do actually form. 
The formation, abundance and properties of halos below 
a cutoff scale has been studied in detail in the context of 
warm dark matt er (WDM) models, both for field halos 
e.g. IBode et~al"ll2001t IKnebe et aLl l 2003| l Yoshida et aT 



has little or no trend with host mass or redshift. The 
figure shows that our SUSY sub-microhalo velocity func- 



2003) and subhalos (e.g. IColm et al.l2000HYoshida et~ 
20031 ). A quantitative comparison with these results lies 



beyond the scope of this work. In this section we illus- 



6 



Early SUSY-CDM substructure 




[km/s] 

Fig. 5. — SUSY sub-microhalo circular velocity function at z = 
86 (solid curve), z = 83, z = 80, z = 78 (long-dashed curves, from 
top to bottom), and z = 75 (short-dashed curve). Volume units are 
comoving. The solid straight line is the universal subhalo velocity 
function of [Reed et al. (2005|), extrapolated from 20 — 300 kms - 1 
down to scales of ms -1 . The dotted lines illustrate deviations by 
factors of four. 



trate one formation path of sub-M c halos and place an 
upper limit on their abundances in order to ensure that 
DM clumps below M c do not provide a significant con- 
tribution to the DM annihilation signal. 

The perturbed cubic grid initial conditions used here 
have a characteristic length scale and preferred direc- 
tions, which is not the case fo r unordered "glass" ini- 
tial co nditions l|White lfl993j) . iGotz fc Sommer-Larsenl 
l)2003|) suggested that grid initial conditions lead to ar- 
tificial fragmentation of WDM filaments and that glass 
initial conditions should be used instead to overcome this 
problem. Even glass WDM simulations show t he for- 
mation of sub-cutoff objects <)Yoshida et al.ll20fKtk how- 
ever, meaning either that glass initial conditions do not 
prevent artificial fragmentation or that the formation of 
sub-cutoff structure is real. Our sub-A/ c halos do not 
form with a preferred separation or along the initial grid 
axis (see Figs. ^andE)), which implies that they are 
not artifacts produced by regular grid noise in the initial 
conditions. 

Figure El shows the formation of a clump of mass 
0.3 x 1CT 6 M Q at z = 86. This halo forms 'top-down' 
through the fragmentation of a fila ment, like the smallest 
WDM halos seen in lKnebe et"aD <|2003l) . Note that the 
fragmentation occurs only after the formation of larger 
halos in the same region, i.e. the sub-M c clump forms 
after its larger neighbors and it therefore has a lower 
peak density and concentration. One can think of the 
sub-Af c clump as a low amplitude initial overdensity on 
some scale larger than M c which fails to grow to a mass 
above M c : by the time the low amplitude perturbation 
goes non-linear its larger and higher amplitude neigh- 
bors have already collapsed and are accreting some of 



the mass which in linear theory would correspond to the 
smaller, low amplitude peak. The small peak collapses 
at a later time in a relatively depleted environment and 
fails to grow above M c . 

Structures above M c = M Q from almost at the 

same time (see section l3~T|) and those below M c appear 
to form even later. This explains the low concentra- 
tion and density contrast of sub-A/ c clumps. Their mea- 
sured abundance depends strongly on the operational 
definition of substructure, i.e. on the subhalo finding 
algorithm used and its parameter values (see Appendix 
for details). The largest abundance of sub-M c clumps 
is found with SKID usin g a density estimate over 64 
particles (see Figure IA9|I . This abundance must be 
taken with caution, however, since many of the smaller 
groups identified in this way do not show up as signifi- 
cant peaks in local density or phase-space density maps. 
We use the s = 64 SKID value only as an upper limit 
to the number sub-M c clumps. It demonstrates that 
the subhalo mass function becomes much shallower than 
dn/dM suhh oc M~^ h below about 0.5 M c (Figure EJ). To- 
gether with their late formation times and low concentra- 
tion this is sufficient to exclude a significant contribution 
to the DM annihilation signal (see section^} from clumps 
below M c . 

3.4. Survival 

In order to properly quantify the survival probabili- 
ties and robustness against tidal stripping of subhalos at 
high and low rcdshift we have traced substructure in the 
entire high resolution region forward in time. We ex- 
amine three time intervals, each covering a comparable 
increase in scale factor. Because of the rapid assembly 
of the SUSY halo we are limited to a relatively small ra- 
tio of final to initial scale factor, a//flj w 4/3, i.e. from 
Zi = 115 to Zf = 86 and from Zj = 101 to Zf = 75 for the 
SUSY halo, and from z; = 0.333 to the present epoch for 
the B9 cluster. Subhalos are identified with 6DFOF (see 
Appendix [SJ and we consider only halos with more than 
320 linked particles at z.;, an order of magnitude larger 
than our usual choice of 32 linked particles at Zf. This 
allows us to identify surviving structures that have lost 
a large fraction of their mass (up to 90 percent) due to 
tidal effects. A halo at z, is counted as a survivor if it 
contributes at least 10 particles to some halo at Zf, unless 
it has merged with the central core of the Zf parent sys- 
tem. Note that even halos for which we find no remnant 
may actually have one at higher numerical resolution. 

We approximate the Lagrangian region of the parti- 
cles that will end up within the virial radius of the par- 
ent at 2/ by a sphere of radius Rc,vh centered on the 
most massive progenitor (MMP) at Zi, and enclosing a 
mass Mhaio(-Zf)- Halos within Rc.vir at Zi are likely to 
end up within r V j r at zj. Table summarizes the sur- 
vival probabilities of well resolved substructure as a large 
host system is forming. Considering all halos in the high 
resolution region (top three rows) we see that early sub- 
structure is significantly more likely to be disrupted, with 
more than 99% of all subhalos surviving from Zi = 0.333 
down to the present epoch, but only 87% surviving from 
Zi = 115 to Zf = 86 and even less (73%) from z, = 101 
to Zf = 75. Tidal destruction is more efficient closer to 
the center of the host, and this causes a further reduc- 
tion in the survival fraction of substructure within two 
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Fig. 6. — Projected DM density maps illustrating the formation of a halo with M^aio = 0.3 X 10 6 Mq, i.e. below the cutoff scale 
of M c = 10 — 6 Mq. Three orthogonal projections of the same cube are shown at each time. The cubes have a size of L = 2.7 comoving 
pc and are centered on those particles which end up in the sub-M c clump at z = 86. The logarithmic color-scale goes from p CI it{z) L to 
10 2 4 p cr it( 2 ) L a t all times. The whole region is relatively overdense: its mean density increases from 6p cr it( z ) at z = 116 to 21p cr ; t (z) at 
z = 86. The three projection directions are aligned with the axes of the initial grid and the initial particle spacing of 0.3 comoving pc is 
indicated by the ruler in the lower-right panel. Note that these structures do not show a preferred direction, alignemcnt with the axes of 
the initial grid, or any regular, grid-scale separation which would indicate artificial fragmentation caused by the grid initial conditions. 



(middle three rows) or one Rc,vir (bottom three rows). 



In Figure [7] we plot the survival fraction versus distance 
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TABLE 2 

Numbers and fractions of surviving structures 



Region 




■-/ 


JVtot 




/sur 


core 




all 


115 


86 


150 


131 


0.873 


3 


1.286 




101 


75 


227 


165 


0.727 


2 


0.954 




0.333 





774 


771 


0.996 


1 


0.956 




115 


86 


111 


94 


0.847 


3 


1.242 




101 


75 


206 


145 


0.704 


2 


0.944 




0.333 





497 


494 


0.994 


1 


0.925 


< Rc, vir 


115 


86 


75 


58 


0.773 


3 


1.033 




101 


75 


121 


75 


0.620 


2 


0.765 




0.333 





261 


258 


0.989 


1 


0.8414 



Note. - Substructure is identified based on phase-space 
density (using 6DFOF) at z, and its survival is checked at ze. 
JVtot and N aur are the total number and number of surviving 
halos, and / sur = N SUT /N tot denotes the fraction of surviving 
substructure. N COIC is the number of halos merging with the 
core of the Zf halo. F Vc max = (f c ,max,i)/("c,max,f ) denotes the 
ratio of the initial to final average maximum circular velocity. 
The first group of three rows is for the entire high resolution 
regions, the second and third groups arc for substructures initially 
within two and one Rc v ; r , respectively, where Rc, vir is the 
radius encompassing M^i (z-r) at Z{. See text for details on the 
6DFOF phase-space halo finder and on how survival is determined. 



to the MMP in units of Rc,vir, ^mmp, for the three time 
epochs analyzed. This drops rapidly for g?mmp ^ Rc,viv 
in the case of SUSY sub-microhalos, below 50% (30%) 
for the epoch ending at Zf = 86 (zf = 75). 

The subhalos that do survive experience significant 
mass loss due to tidal stripping. We quantify this mass 
loss by calculating the ratio of the average of jj c max of 
all halos at the beginning to that at the end of each 
epoch, F„ cmax = (w c ,max,i)/(wc,max,f) (Figure EJ. For ha- 
los beyond Rc, vir ^c.max grows with time in the SUSY 
simulation, whereas for B9 the ratio remains about con- 
stant. The increase in average w Cjma x is especially strong 
in the redshift interval 115 — 86, F Vcmtac w 1.5, and this 
again illustrates the rapid and late (in terms of a//aj) 
formation and growth of these structures. SUSY sub- 
microhalos within i?£, v ir suffer more mass loss than those 
beyond R c , vi r, i^ , m J(< Rc,™) = 0.68^ max (> i? £ , vir ). 
For the cluster substructure the difference between outer 
and inner regions is less pronounced, indicating that it 
is more resistant against tidal stripping. This robustness 
translates into a higher survival probability, as shown 
above. 

4. DARK MATTER ANNIHILATION SIGNAL 

The differences in substructure contrast described in 
the previous sections have ramifications for the expected 
7-ray signal from DM self-annihilation. Indirect detec- 
tion experiments are hoping to detect both the cosmolog- 
ical extragalactic background as well as the local signal 
from the Milky Way's center and from DM substruc- 
ture in its halo. Since the annihilation rate is propor- 
tional to p|) M , the predicted flux depends sensitively on 
the dumpiness of the DM distribution. Small-scale sub- 
structure is expected t o significantly increase the signal 
from individual hal os llC alcaneo-Roldan fc Moorell2000t 
IStoehr et all 12003 IColafrancesco et alJ I2005ir~as well 
as the the diffuse extragalactic background (lUllio et alJ 
l200l iTavlor fc Silkll200a lAndo fc Komatsull2006|) . An- 
nihilation luminosity estimates based on DM densities 




1 2 



Fig. 7. — Fraction of surviving structures (solid circles) and 
ratios of initial to final average circular velocity i> c ,max of surviving 
substructure (open squares) versus initial separation from the most 
massive progenitor in units of Rc v ir- Halos were traced from z, = 
115 down to Zf = 86 (solid lines), from z, = 101 to zj = 75 
(dashed), and from Zj = 0.333 down to the present epoch (dash- 
dotted). 



measured in collisionless N-body simulations of a Milky 
Way-scale halo suggest that substructure may boost the 
signal by only a small factor, ~ 70%, compared to a 
smooth spherical NFW halo (Stoehr et al. 2003, here- 
after S03). In this work we follow S03's analysis and 
compare the contribution of substructure in the SUSY 
host at z = 75 and the B9 galaxy cluster at z = 0. 

The annihilation luminosity is proportional to the ra- 
tio of the velocity weighted self-annihilation cross section 
and the mass of the DM particle, (av)/mx- Following 
S03, we disregard this (highly uncertain) particle physics 
term and concentrate instead on the astrophysical com- 
ponent, defining the annihilation signal of the i th halo to 
be 

Si = / &u dV i = Y] Pj m j, ( 2 ) 
Jv < MP,} 

where pj and mj are the density and mass of the j th par- 
ticle, and {Pi} is the set of all particles belonging to halo 
i. Whereas S03 estimated the pj with the parameter- free 
Voronoi tessellation technique, we use the SPH-kcrncl es- 
timated densities from SMOOTH over 64 nearest neigh- 
bors. 

Like S03 we find that the component from the 
spherical-average density profile contributes only a frac- 
tion of the total signal, 49% in the SUSY halo and only 
35% in the B9 halo. The contribution to the signal from 
all resolved bound subhalos is 39% in B9 and only 10% in 
the SUSY host. The remaining difference is due to both 
unbound streams and departures from spherical symme- 
try 4 . The absolute contribution of the resolved substruc- 

4 These estimates depend only weakly on the smoothing scale 
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ture component depends on numerical resolution, since 
higher resolution simulations would be able to resolve 
highe r densities in the subhalo centers (|Kazantzidis et all 
|2004|) . The bound substructure contribution to the total 
signal could be larger and the values given here should be 
taken as lower limits. Both the SUSY and B9 halos, how- 
ever, have approximately the same number of particles 
within their virial radii, allowing us to directly compare 
the contribution from resolved substructure between the 
two. We conclude that it is indeed the lower density con- 
trast in the early SUSY host that leads to a significantly 
smaller contribution from bound sub-microhalos to the 
total annihilation signal. 

used for the SPH density estimates: using 1024 instead of 64 
nearest neighbors gives relative contributions of spherical smooth 
signal, bound substructure, and remaining density fluctuations of 
(57%, 9%, 34%) instead of (49%, 10%, 40%) for the z = 75 SUSY 
halo. 
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Fig. 8. — The ratio of the annihilation signal of individual sub- 
halos to the total signal from all dark matter particles in the host 
halo versus the subhalo mass ratio M^^/M^^. Circles denote 
the signal calculated directly from the subhalo particle densities, 
triangles are used for an estimate of the signal assuming an NFW 
profile and a concentration c PS 2r v j r /r Vc max . The top panel shows 
the B9 halo at z = 0, the bottom one the SUSY halo at z = 75. 
Subhalos with the same signal-to-mass ratio as the host would lie 
on the thick solid line. The dashed lines show the best-fit linear 
relation between subhalo signal and normalized mass (lower line: 
direct calculation, upper line: NFW estimate). Halos to the right 
of the thin vertical line have more than 200 particles. 



In Figure |S1 we present the mass dependence of the 
annihilation signal for individual subhalos. To ease a di- 
rect comparison we normalize the subhalo mass by the 
total halo mass and the subhalo signal by the total sig- 
nal from all halo DM particles. Since the signal strength 
may depend sensitively on the internal density structure 
of the subhalo, we only include halos with more than 200 
particles in our analysis. If substructure had the same 
signal-to-mass ratio S/M as their host they would lie on 
the thick solid line in the figure. We find instead that 
the B9 subhalos tend to have a larger S/M (i.e. they lie 
above the solid line, in qualitative agreement with S03), 
but that the majority of the SUSY sub-microhalos have 
lower S/M than their host. Furthermore we find that 
S/M is not independent of mass: in both our simula- 
tions less massive substructures have a lower S/M than 
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more massive ones. For B9 we find S/M oc M^L, while 
for SUSY S/M oc M° u £ h . 



ubh' 

Given a substructure abun- 



dance of dN/dM SX ibh oc M~?jj with a < 2 (see Section 
I3.2f> . this suggests that the total signal may be domi- 
nated by the highest mass subhalos. One has to keep in 
mind, however, that the peak densities in an N-body halo 
are always limited by numerical resolution. The local 
dens ities are too low both in the center of the host sys- 
tem Ipiemand et al.l[2004fl and even more so in the inner 
regions of subhalos i Kazantzidis et al"1l2004|) . It is thus 
likely that estimates of the total annihilation signal based 
directly on densities measured in N-body halos underes- 
timate the importance of the substructure contribution . 
Analytic calculations (|Calcaneo-Roldan fc Moorell2000j) 
tend to find larger substructure boos t factors and a signal 
that is dominated by small subhalos (jColafrancesco et alJ 

We have attempted to correct for the artificial heating 
of the densest, inner regions by recalculating the annihi- 
lation signal for every subhalo assuming an NFW profile 
with a concentration given by c w 2.16r v i r /r Vc max . In 
this case the normalized annihilation signal is equal to 



S/St, 



1.22 



G 2 r v 



/St, 



(3) 



where G is the gravitational constant, and w c>max and 
r Vc max are determined using SKID with spherical binning 

of the bound subhalo particles. Expressing Si in terms 
of i> c ,max an d r Vc max reduces the dependence of the esti- 
mate on the strongly resolution dependent inner density 
profile. These corrected subhalo signals are plotted as 
triangles in Figure |SJ As expected the NFW estimate 
enhances the signals, and the correction is about 30% 
in both the SUSY and the B9 halo. With this boost 
the relative contributions of spherical smooth signal, 
bound substructure, and remaining density fluctuations 
are (47%, 14%, 39%) for the SUSY and (30%, 48%, 22%) 
for the B9 halo. Note that the substructure signal in the 
B9 run is still significantly larger than in the SUSY run. 
The mass dependence of S/M appears to be preserved, 
but we caution that even our NFW estimates are not bias 
free: the finite resolution in our sim ulations may well 
lead to an underestimate of i> c ,max (c£ IKazantzidis et all 
|2004|) . especially for low-mass subhalos. 

Our results suggest that in a real sub-solar mass SUSY 
halos the signal from subhalos and other density fluctua- 
tions is at least as big as the signal from the spherically- 
averaged component. Estimates o f the cosmologi- 
cal, extragalactic 7-ray background |uTlio et alJ 120021: 
iTavlor fc Silkll2002t lAndo fc Komatsul l2006j) should be 
revised to take i nto ac count both microhalos in the field 
;U1l2005D and sub-microhalos within larger 



IjDiemanc 

(sub)-solar mass systems like the one studied in this 
work. We plan to address remaining open questions 
(i.e. the total substructure boost factor and what sub- 
halo mass scale dominates the signal) in future work us- 
ing higher resolution simulations, extensive convergence 
tests, and comparisons with various analytical estimates. 

5. CONCLUSIONS 

We have simulated the formation of a 0.014 M Q halo 
at z = 75 in a scenario where CDM consists of 100 GeV 
supcrsymmctric particles. Our main results are: 



• Early SUSY-CDM structures form almost simul- 
taneously over a wide range of mass scales. Halo 
masses grow very rapidly, and filaments are very 
extended compared to the size of the largest halos. 

• The contrast between substructure and back- 
ground is lower in early SUSY halos than in galaxy 
cluster hosts, both in density and in phase-space 
density. 

• At z = a simple, fast phase-space friends-of- 
friends subhalo finder (6DFOF) gives robust re- 
sults. Velocity functions and subhalos number den- 
sity profiles are identical to those obtained with 
SKID. Unfortunately on SUSY scales both meth- 
ods have problems due to the reduced contrast. 
Results from both methods depend quite strongly 
on the input parameters, but can be brought into 
agreement with a suitable choice of parameters. 

• In relative units both the sub-microhalo mass and 
circular velocity functions of the 0.014 M Q SUSY 
host are in the range found for z = galaxy clus- 
ters. 

• In absolute units both the circular velocity func- 
tions of the 0.014 Mq host halo lie within the scat- 
ter of the univ ersal subhal o circu lar velocity func- 
tion found by iReed et all l)2005H in low redshift 
dwarf, galaxy, and cluster halos. 

• As the universe expands by a factor of 1.3, about 20 
to 40 percent of well-resolved early substructure is 
destroyed, compared to only about 1 percent in the 
low-redshift cluster. The reductions in the peak cir- 
cular velocities among the survivors are also larger 
for the early micro-subhalos. 

• The contribution to the DM annihilation sig- 
nal from the bound sub-microhalos resolved in 
our simulations is about four times smaller than 
the signal from the spherically-averaged compo- 
nent in the SUSY host. In a cluster simulated 
at similar numerical resolution (B9), the bound 
substructure contribution exceeds the spherically- 
symmetric part. Due to numerical limitations these 
are lower limits. The reduced contrast in early- 
forming sub-microhalos results in a lower substruc- 
ture boost factor compared to galaxy or cluster ha- 
los. 
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APPENDIX 

IDENTIFICATION OF EARLY SUSY-CDM MICRO-SUBHALOS 

The simultaneous nature of structure formation near the SUSY cutoff scale at early times leads to substructure 
with a reduced contrast both in real and phase-space density. Conventional subhalo finders use positional information 
in the group finding stage and velocities only after the fact to decide which particles are gravitationally bound to a 
subhalo. In this appendix we discuss how a conventional subhalo finder (SKID) performs in these challenging low 
contrast systems and present a simple and fast new method of finding (sub-)structures in phase-space. 

SKID: self bound over-densities 

SKID (publicly available at: www-hpcc.astro.washington.edu/tools/skid.html) is described in detail in lStadej l)200lF) . 
In a dark matter only simulation the group-finding algorithm goes through the following steps. First, the lo- 
cal densities are calculated using an compact SPH kernel whose scale is adjusted for each particle to includes 
a fixed number V of nearest neighbors (the same method is used in the program SMOOTH available at www- 
hpcc.astro.washington.edu/tools/smooth.html). Second, particles are moved along the density gradients until they 
move less than a user specified length e, i.e. they are trapped in a local potential well of scale e. When all particles 
satisfy this criterion, they are linked together using a linking length equal to e. Finally the gravitationally unbound 
particles are removed iteratively (one by one) from all the groups. 

SKID results depend on two parameters: the linking length e and the number of nearest neighbors s. A small e tends 
to underestimate the mass of larger subhalos while a large e links nearby subhalos together. In current high resolution 
z = CDM clusters (roughly ten million particles within the virial radius) setting e = 0.004 r v ; r successfully identifies 
most resolved structures. Combining subhalo catalogs derived wit h a range of values for e results in small c orrections 
at the low mass and high mass end of the subhalo mass function ijGhigna et al.ll2000l iDiemand et al.ll2004|) . For the 
z — 86 SUSY halo we find a similar weak e-dependence of the SKID results. 

The number of nearest neighbors s is usually set to s = 64 and in a z = CDM cluster the SKID subhalo mass 
and velocity functions are practically independent of s, with the exception that some subhalos with less than about 
0.25s bound particles are lost when a larger smoothing scale is used (see Figures lA~9l and |A1 111. Unfortunately SKID's 
results are much more parameter dependent for the SUSY halo: the right panel of Figure lA9l shows that the SUSY 
subhalo mass function depends on s over the entire subhalo mass range. Many of the small subhalos found with s = 64 
and s — 256 appear spurious, i.e. visual inspection fails to identify a clear corresponding peak in the mass (or phase 
space) density projections. Only at the high mass end there is hint of convergence when using s > 1000. Tests in the 
z = cluster suggest that using large s leads to an underestimates of the subhalo abundance below about 0.25s bound 
particles. We compromise and use s = 1024 throughout the paper, but caution the reader that the resulting group 
catalogs are probably incomplete below about 256 bound particles. 




Fig. A9. — SKID subhalo mass functions for the z = galaxy cluster (left) and the z = 86 SUSY CDM halo (right) using different 
numbers of neighbors s for the SPH density calculation. Plotted are s = 64 (solid), 256 (dotted), 1024 (dashed), and 2048 (dash-dotted, 
SUSY halo only). The first bin always begins at 10 bound particles. Bars show the Poisson errors. For the massive halo, the mass functions 
are independent of s above about 0.25s. Due to the smaller density contrast the entire mass function depends on s for the microhalo. 
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Finding structures in phase-space: 6DF0F 

Subhalos have smaller internal velocity dispersions a than the hot uniform density background. Figure[21show maps 
of p/cr 3 (both p and a are calculated using an SPH kernel over 32 nearest neighbors), which is a simple and fast way 
to estimate the local phase-space density. The phase-space density contrast is much larger than the real-space density 
contrast. Even small subhalos, which due to limited numerical resolution don't reach peak mass densities as high as 
their host halo, stand out more clearly in phase-space density because internal velocity dispersions are smaller for less 
massive halos. In order to take advantage of the increased contrast in phase-space we have developed a new, simple, 
and fast (sub-) halo finder. 

The algorithm is based on the friends-of-friends algorithm (FOF) ({Davis fc Diorgov ski 1985), in which particles are 
linked together if they are separated by less than some linking length b x dx, where b < 1.0 and dx is the mean particle 
separation. We simply extend the standard FOF algorithm to include also a proximity condition in velocity space. 
Two particles are linked if 

(si - x- 2 ) 2 /(b x dx) 2 + («! - v 2 ) 2 /b 2 v < 1 , (Al) 

where b is the standard FOF linking length in units of the mean particle separation dx and b 2 is the velocity space 
linking length. For b 2 — > oo one recovers the usual real-space FOF. After the linking stage all groups with at least 
N min are kept and circular velocity profiles are calculated around the center of mass of each group using 100 linear 
spherical bins. Since 6DFOF finds only the peaks in phase-space density we set the radius of the outermost radial bin 
to 5 times the half mass radius of the linked group. The steps described above were implemented in parallel within 
PKDGRAV and they run at with vanishing computational expense compared to the gravity calculations. We omit the 
computationally expensive step of removing unbound particles and hence we don't actually determine subhalo masses. 
Instead we estimate v c ,max as a proxy for mass by fitting the circular velocity profiles using the circular velocity 
corresponding to an NFW subhalo embedded in a constant density background. The three free parameters in these 
fits are u c ,max and r Vc max of the subhalo and the background density. We minimize the absolute velocity differences. 

Additionally we only allow fits where the subhalo mass enclosed within r Vc max is larger than the background mass. 
This constraint prevents a few cases (about 2% of the groups) where small fluctuations in the background density at 
large r adii a re fitted with an extended, diffuse NFW halo which lies well below the background density at most radii. 
Figure lAlOl shows some example circular velocity profiles around 6DFOF-structures found within the virial radius of 
the z = 86 SUSY halo. Typical z=0 subhalos are much denser than the local background and show a clear peak in the 
total circular velocity profile. Early subhalos, on the other hand, often lack a clearly distinguished pea k and appear 
as excess circular velocity above the linear circular velocity from a constant density background (Figure lAlOfl . 

After some experimentation we found that with b = 0.04, b v = 295 kms -1 and iV m i n = 32 the SKID subhalo velocity 
function and radial distribution for the z = cluster B9 are almost exactly reproduced. They differ significantly 
only below about 0.025t ' n.Tna->t.host. which lies below the complete ness limit of about 0.05?; Cimax .host in such a simulation 
fe.g. lGhigna"eTld1l2^0rllDiemand et alll200l iReed et alJl2005|) . The 6DFOF results are not 'sensitive to small (15%) 
changes in both b and b v . The adopted velocity linking of b v = 295 kms -1 is about four times smaller than the mean 
ID velocity dispersion within the virial radius. This suppresses the random links between background particles even 
in dense regions were each sphere of radius 0.04 x dx contains many particles. Only a prolate structure enclosed in 
the inner 15 kpc (about 1.5 percent of the virial radius) is linked together. In this paper we don't count this central 
structure as a subhalo 5 . 

In the SUSY halo at z = 86 we chose a somewhat larger space linking length due to the smaller peak overdensities 
b = 0.057 and rescale the velocity linking lengths to b v = 4.15 ms _1 (four times smaller than the mean ID velocity 
dispersion within the virial radius) and we keep N m i n = 32. Unfortunately just as with SKID also the 6DFOF subhalo 
lists at z = 86 depend strongly on the choice of group-finding parameters used. A 15 percent increase (or decrease) in 
both b and b v yields a very different velocity function, whereas at z = our 6DFOF results are unaffected by the same 
changes in b and b v (Figure lAlljl . Reducing the 6DFOF linking lengths by 15 percent (a factor of 2.3 in phase-space 
density threshold) results in 40 percent less subhalos (459 vs. 779) within the virial radius and the resulting subhalo 
velocity function is about a factor of five lower. Increasing the linking lengths by 15 percent results in an increase of 
68 percent in the number of subhalos (1312 vs. 779), but many of them have maximum circular velocities well below 

0. 3 mr 1 , and show no clear subhalo signal in their total circular velocity profile v c (r) = \J v t su b + v c bg- Visual 

inspection of these small v Cj max , sub structures also shows no convincing peak in either real or phase-space density. 
Only for the highest w c . max, sub structures do the longer linking length results agree with those using our fiducial 
values of b = 0.057 and b v = 4.15 ms _1 . This strong parameter dependence of 6DFOF results at z = 86 is again a 
consequence of the similar formation times of small mass, high redshift subhalos and hosts which reduces the contrast 
in phase-space density. 

When tracing structures over time we set bx dx and b v to constant values in proper coordinates. This corresponds to 
using the same proper minimum phase-space density to define structures at different epochs. This choice is motivated 
by the fact that the high (phase-space) density peaks we aim to identify have long decoupled from the Hubble expansion, 

1. e. they correspond to the inner cores of virialized systems. Employing a constant comoving linking length would 
result in a decreasing phase-space density threshold. The values we use are bxdx = 0.04 x 0.111 kpex 1/(1 + 0) = 4.444 
proper kpc, b v = 295 kms^ 1 for the B9 cluster and b x dx = 0.057 x 0.3 pc x 1/(1 + 86.3) = 0.197 proper Mpc and 



5 SKID produces a similar central object which is always discarded in this work as well. 
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Fig. A10. — Example circular velocity profiles (solid) of structures found with 6DFOF within the virial radius of the z = 86 SUSY CDM 
halo. Best fit total circular velocity profiles v c (r) (dashed- dotted) contain a contribution from an NFW subhalo v c su b (dashed) and from a 

homogeneous background density v c -qq (dotted), i.e. v c (r) = sub + wj? BG ■ Upper panels show small subhalos below the completeness 

limit of about 0.05v c maXj host = 1 ms _1 and the lower panels show larger subhalos. From left to right subhalos lie closer to the host center, 
i.e. the background density increases. In two cases (middle panels of the top and bottom rows) we also show the best fits obtained without 
the constraint that the subhalo mass must be larger than the background within r Vc max (long- dashed). Velocities and radii are given in 
physical units. 



b v = 4.15 m4 for the SUSY halo. 
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finder was run with b = 0.04, b v = 295 kms -1 and iV m ; n = 32 (thick solid) and 15% larger (upper solid) and 15% smaller (lower solid) 
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